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Abstract 

We study geodesies on the parameter manifold, for systems exhibiting sec- 
ond order classical and quantum phase transitions. The coupled non-linear 
geodesic equations are solved numerically for a variety of models which 
show such phase transitions, in the thermodynamic limit. It is established 
that both in the classical as well as in the quantum case, geodesies are con- 
fined to a single phase, and exhibit turning behavior near critical points. 
Our results are indicative of a geometric universality in widely different 
physical systems. 
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1 Introduction 



Information theoretic studies of phase transitions are, by now, well established. 
The underlying idea here is geometric in nature, and rests on the definition of a 
Riemannian metric tensor on the space of parameters (called parameter manifold) 
of a system. Depending on whether the interactions of the system are classical 
or quantum in nature, this metric might be induced from the equilibrium ther- 
modynamic state space [1] (for a review, see pj), or from the natural Hilbert 
space structure of quantum states [3]. For the former, the parameter manifold 
consists of thermodynamic control parameters such as the pressure, volume and 
temperature, while for the latter, this might be thought of as the manifold of 
coupling constants appearing in the Hamiltonian. 

Given such a metric tensor, the parameter manifold can have very different 
properties depending on whether the system undergoes a second order classical 
or a quantum phase transition (CPT or QPT). Whereas the hallmark of a CPT is 
that the scalar curvature arising out of the metric diverges at a second order phase 
transition (and everywhere on the spinodal curve), this is not the case for second 
order QPTs where the curvature can remain regular [4|- It is also known that 
whereas some components of the metric tensor vanish at a second order CPT, 
as these are related to inverses of thermodynamic response coefficients [2], for 
QPTs, the situation is reversed, and some of the components of the metric tensor 
diverge at such a transition, as follows from first order perturbation theory |1] 
(although this may not be true in some special cases, see [5]). 

Although a lot of attention has been paid to the behavior of the metric tensor 
and its associated scalar curvature in the context of phase transitions, much less is 
known about geodesies, i.e paths that minimize the distance between two points 
on the parameter manifold. In any geometric setup, the behavior of geodesies 
is an important object to study. Some studies on geodesies have appeared in 
the context of CPTs [6], and QPTs (specifically, for adiabatic quantum compu- 
tation) [7], [8], in special cases. The purpose of this paper is to complement 
and generalize these results, and to obtain and analyze general solutions to the 
geodesic equations for some model systems that exhibit second order phase tran- 
sitions. 

Here, we study four models in the thermodynamic limit : the Van der Waals 
(VdW) model for fiuids, the Curie Weiss (CW) mean field model of ferromag- 
netism, the infinite Ising ferromagnetic chain - all of which exhibit CPTs at 
finite temperature, and the transverse field XY model that exhibits a QPT at 
zero temperature. For all these models, the full set of coupled non-linear geodesic 
equations in the information geometric context are set up and solved numerically, 
with appropriate initial conditions. To the best of our knowledge, such an analy- 
sis has not been performed before. Our treatment is completely general in nature, 
and differs significantly from the methods used in [6] , [8] where the focus was on 
obtaining specific geodesies between two given points in the parameter manifold. 
Interestingly, we find that in all the examples that we consider, geodesies exhibit 
a "turning point" close to criticality, and are "confined" to a single phase, thus 
indicating a geometric universality in apparently unrelated physical phenomena. 

This paper is organized as follows. In the next section, we first briefiy recall 
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some basic facts about information geometry and geodesies. We then proceed 
to analyze the VdW, the CW and the infinite Ising ferromagnet, as illustrations 
of CPTs. For the Ising ferromagnet, our analysis of information geometry is 
novel, and has not appeared in the literature before. In section 3, we analyze the 
geodesic structure of QPTs via the transverse field XY spin chain. We end in 
section 4 with our discussions and directions for future study. 



2 Information Geometry, Geodesies, and Clas- 
sical Phase Transitions 

In the context of equilibrium thermodynamics of classical systems, the formula- 
tion of information geometry is mainly due to the work of Ruppeiner |2]. The 
main idea here is to consider the positive definite Riemannian metric arising out 
of the Hessian of the entropy density s, and given by a line element 

Here, x'^,// = 1,2, denotes the internal energy and the particle number per unit 
volume, and are co-ordinates on the parameter manifold in the "entropy repre- 
sentation." ks is the Boltzmann's constant, which we will set to unity in what 
follows. The line element of eq.([T]) introduces the concept of a distance in the 
space of equilibrium thermodynamic states via fluctuation theory, i.e, the larger 
is this distance between two given states, the smaller is the probability that these 
are related by a thermal fluctuation. Various representations (related to each 
other by Legendre transforms) can be used for this geometric construction (a full 
list can be found in [2]), and a particularly useful diagonal form of the metric for 
single component fluids and magnetic systems is 

where T is the temperature, p the number density, and fi = > / being 

the Helmholtz free energy per unit volume. For magnetic systems, we need to 
consider thermodynamic quantities per unit spin, with the magnetization per 
unit spin m replacing p. 

On the other hand, information geometry in quantum mechanical systems, 
first studied by Provost and Vallee [3j, is defined by taking two infinitesimally 
separated quantum states and constructing the quantity 

\ip {x + dx) — ijj (x) p = {d^il)\dyil))dx'^dx^ = a^^^ydx^dx^ (3) 

where x^ (collectively denoted as x in the l.h.s of eq.dH])) denotes the parameters 
on which the wave function ip depends on, and (9^ is a derivative with respect to 
x^. From the a^i, (which are not gauge invariant), a meaningful gauge- invariant 
metric tensor can be defined as [3] 

QtJiu = o^^lu - l^^ll^v] {3^ = {x)\d^i) {x)) (4) 
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Here, gfj,,y is the metric induced from the natural structure of the Hilbert space 
of quantum states. The metrics in eqs.Q and (jl]) can be used to predict second 
order phase transitions in both CPTs ^ and QPTs ^4j. We also record here the 
expression for the scalar curvature arising out of the metric in the special case 
when the metric is diagonal (with g = det g^^iy) : 
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(5) 



Given the information geometry of classical or quantum systems, we wish to 
study geodesies in the same. Let us briefly recall a few elementary facts about 
geodesies. For a manifold endowed with a metric with components (^^j,, a geodesic 
is a path that extremizes the proper distance (or line element, whose infinitesimal 
form is given by dr'^ = g^^dx^dx'^). This can be cast as a variational problem, to 

determine the extrema of the integral ^ g^j^^xi^x^dX where the dot denotes a 
derivative with respect to A, which is an affine parameter, parametrizing the curve 
joining two points denoted 1 and 2. Calculus of variations can then be applied 
with the result that geodesic curves are solutions to the differential equations 

*- + a«^^o. with r;,49-c(%. + |j_|^) (6) 



The above equation can also be obtained by writing a "Lagrangian 

1 
2 



£ = - {g^.x^x'') (7) 



and using the (derivatives of the) Euler-Lagrange equations that follow. This 
method often provides valuable insights into the symmetries of the system. We 
will be interested in studying the solutions of eq.(l6]) in the context of CPTs and 
QPTs. It is well known that a natural affine parameter for geodesic curves is 
A = r, and thus it is useful to consider the normalized vector x^ = dx^/dr such 
that x^x^ = g^yX^x^ = 1. 

Eq.dn]), in general, gives rise to a set of coupled non-linear differential equa- 
tions, which might be difficult to solve analytically. We will mostly concentrate on 
numerical solutions with appropriate boundary conditions. Note that in terms of 
the normalized vector x^^, we need to specify three boundary conditions in order 
to solve eq.([6]), with the fourth one being fixed by the normalization condition. 
Namely, we choose a "starting point," i.e, an initial value of x^, and any one 
component of x'^. The second component of the derivative is then determined 
from the fact that x^ is normalized. With these boundary conditions, we deter- 
mine the most general solutions to eq.(l6]), and study geodesies near criticality. 
This is done by solving for x^ in terms of the affine parameter r, and tracing 
out the geodesic near the critical point, by parametrically plotting the resulting 
solution, under variation of r. Let us now illustrate the above discussion with 
the example of the Van der Waals fluid and the Curie- Weiss ferromagnet. 



^We will also keep in mind that geodesic paths are not unique : an elementary example is 
that of a 2-sphere, where there are an infinite number of geodesies, i.e great circles, between 
two anti-podal points. 



3 




Figure 1: Numerical solution for geodesies of the VdW equation of state close to 
criticality, in the {pr,Tr) plane. The dashed blue, dot-dashed red and solid green 
curves correspond to the boundary conditions (Tr,Pr,Pr) = (1.001,1.009,-1.2), 
(1.001, 1.007, -0.92), and (1.0007, 1.011, -2.2) respectively. The geodesies turn 
back from the critical point, {pr,Tr) = (1, 1). 



2.1 The Van der Waals and the Curie- Weiss Models 

Information geometry of the Van der Waals fluid is well established, see e.g 
We start from the Helmholtz free energy per unit volume 



/vdw = -pTln + pc^Tln (^^^ - pTln (1 - bp) - ap^ 



(8) 



where e denotes the exponential function, is the specific heat at constant vol- 
ume, p and T are the number density per molecule of fluid and the temperature, 
and a, b are the coefficients arising in the VdW equation of state. It is conve- 
nient to work with the reduced VdW equation of state, and we can substitute 
a = 9Tc/Spc, b = l/3pc, where pc and Tc denote the critical values of the density 
and the temperature, respectively. Further, the reduced density and temperature 
are defined by pr = p/pc, % = T/Tc. We will set pc = Tc = I to simplify the 
algebraic details, and also choose c„ = 3/2, the ideal gas value. The information 
metric (in terms of the co-ordinates and pr) is then given, from eq.(l2]), by 

3ft _ 9 [47-,, - (p, - 3)'] 

9..- 4p J, (p. - 3)= 

Since we are interested in geodesies close to criticality (for a recent related 
discussion, see |10]), we now expand the metric upto first order about the critical 
point, {Tr,pr) = (1,1) (remember we have set {Tc, pc) = (1,1))- The metric 
components are then given by the simple expressions 

9TT = ^Pr(l-TX g;^ = l{%-l) (10) 
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Figure 2: Various numerical solution for geodesies of the VdW equation of state 
in the {pr,Tr) plane. The dashed blue, dot-dashed black and solid red lines are 
geodesies that begin from the gas, liquid and supercritical phases, respectively. 
The spinodal curve is shown in dotted green. 



where the superscript c in eq. fllOp signifies that these expressions are valid close 
to criticality. The geodesic equations of eq.(l6]) turn out to be 



2Tr - 3 Pr Apr (2T^ - 3) 

I?(27^_ (11) 

- 1 3 (T^ - 1) ^ ' 

We now numerically solve eq.f lTT|) with three boundary conditions : (Tr, Pr, Pr) = 
(1.001, 1.009, -1.2), (1.001, 1.007, -0.92), and (1.0007, 1.011, -2.2). For all the 
three cases, we solve eq.f fTT]) for values of the affine parameter between zero and 
0.0025. The solution for Tr and pr are then parametrically plotted by varying 
the affine parameter. The result is shown in fig. (II]) in the {pr,Tr) plane, where 
the dashed blue, dot-dashed red and solid green curves correspond to the three 
boundary conditions described above, respectively. We see that the geodesic 
curves "turn back" from the critical point. As we will see, this is a generic 
feature for all second order phase transitions studied in this paper. 

For the sake of completeness, we mention here that the analysis of geodesies 
using the full VdW metric of eq.(l9]) is similar, although the geodesic equations 
are more complicated and we omit them for brevity. After extensive numerical 
analysis, our conclusion here is that a geodesic starting in the liquid (p^ > 1, < 
1) or gas {pr < l,Tr < 1) phase does not reach the other phase. They either 
terminate at the spinodal line or continue to the supercritical region. Also, close 
to the critical point, geodesies show the turn-around behavior as depicted in 
fig. (II]). We also find that geodesies do not show any special behavior at the 
binodal lines, i.e at the location of the first order phase transitions, which is 



^The value of Tr is fixed from the normalization condition as alluded to before. 
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expected, because the metric and the scalar curvature are both regular here. 
These results are summarized in fig. ([2]), where we have shown several numerical 
solutions to the geodesic equations for the VdW equation of state. The dotted 
green curve is the spinodal curve. The dashed blue curves on the left and the 
dot-dashed black curves on the right are geodesies that start from the gas and 
liquid phases respectively, and continue into the supercritical region. The solid 
red curves are geodesies in the super-critical region (T^ > 1), and show turning 
behaviour similar to that depicted in fig. ([I]). 

We now move on to study geodesies in the classical mean-field Curie- Weiss 
ferromagnetic model in the thermodynamic limit. Information geometry of this 
model has been studied extensively in [11], and we simply state the result that 
in the (T, m) representation, the line element of eq. (j2]) is given by 

. c.,r + (12) 

J ^ I — 1 

Here, T is the temperature, its critical value, m is the magnetization per unit 
spin, and Cl{T) is a "lattice specific heat" introduced in pjj, that corresponds to 
the mechanical energy of the lattice. As was shown in [12], information geometry 
in the CW model cannot be defined without introducing this term ad hoc in the 
theory. In [Tlj, it was shown that the line element in eq.f ll2p correctly reproduces 
all the known features of the CW model, including the first order phase transi- 
tions. We will study the model close to criticality, and approximate the metric 
close to m = as 

c _ C'l(T) c _ 1 -^c /^ox 

9tT — 7^ ; 9mm ~ ^ 

where again the superscript c denotes that we are close to criticality. To analyze 
the geodesic equations here, we note that a crucial simplification is possible, 
since none of the metric components in eq. f|T3|) depend on the magnetization. 
This implies that the Lagrangian of eq.Q is independent of m, and hence the 
Euler-Lagrange equation that follows from it implies that m = K/g'f^^ where K 
is a constant. Then from the normalization condition g^jT"^ + g'^^m^ = 1, it 
follows that 

^ K^\_ T'[T{l-K^)-n ^^^^ 



9'tt\ 9'mm) CLiT)iT-T, 

It is enough for us to consider the region T > Tc, for which eq.f ll4p implies that 
positivity of the right hand side imposes the restriction T > Tc/(1 — K'^), with 
K'^ < 1. This means that a geodesic in the region T > always remains in that 
region and cannot cross-over into the region T < T^. A pathology arises for the 
case K = 0, for which eq. (fT^ implies that such a restriction is not implied, since 
is always a positive number for = 0, or m = constant. We have checked 
this by explicitly solving the geodesic equations, which in this case are given by 

T2 ( TCr - 2Cr 
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Numerical analyses (after choosing an appropriate regular functional form for 
Cl{T), such as a power series) reveals that geodesies with m = constant lines 
(these are indeed geodesies as they satisfy the second equation in eq.llTS!)) cross 
over inside the spinodal region. This is probably a mathematical artifact and we 
do not have a physical explanation for this. Apart from these constant m lines, 
the behavior of geodesies close to the critical point is, as expected, qualitatively 
similar to that of the VdW fluid, and graphically, they resemble the ones shown in 
fig.([T]). We also find that the behavior of geodesies with the full CW metric (away 
from criticality) is qualitatively similar to those of the VdW model. Specifically, 
geodesies in the phase m > do not reach the phase m < 0, and vice versa. 



2.2 The Infinite Ising Ferromagnet 

We now study geodesies in the infinite-range ferromagnetic Ising model with a 
transverse magnetic field. This model was originally studied in [13], where it 
was shown that in the thermodynamic limit, it can be described by the classi- 
cal dynamics of a single large spin. The information geometric aspects of this 
model has not been studied so far, and we begin with a discussion on this. The 
Hamiltonian for this model is given by [13], [H] 



= E ^^^1 - ^ E = (^tot)' - hs:,, (16) 

i<j 



where the second equality follows from defining the total spin, S^^^ = 
^tot = l^i (and neglecting a constant term). We will set J = 1 in what follows. 
In a mean- field approach, where the average magnetization m = Yli < *S'f > /N, 
the Hamiltonian for a single spin reduces to ifj\p = —mS^^^ — hS^^^. This is an 
effective two-state model whose partition function can be shown to be given by 
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2Co* I I (17) 

To understand the geometric aspects of this model, we write the Gibbs free en- 
ergy for the single spin, G = —TlnZ and effect a Legendre transform to obtain 
the Helmholtz free energy F = G + where m should be thought of as the 

applied magnetic field, i.e an intensive thermodynamic variable. The factor of 
1/2 in the Legendre transform might look strange, but note that this enforces 
the magnetization OF /dm = (via the relation m = —dG/dm), i.e defines the 
boundary between the ferromagnetic and paramagnetic regions. In (T, m) coordi- 
nates, using the expression for the Helmholtz free energy, the metric components 
are given from eq.([2]) by 

9tt = {h^ + m^) Sech^a 

1 1 2 (mV/i^ + m2 + h^TSmh{2a)) ^ ^ 
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Figure 3: Numerical solution for geodesies of the infinite range Ising ferromagnet 
at h = 0.2, near the critical point (T = 0.236, m = 0). All geodesies are chosen 
to pass through the point (T, m) = (0.3,0.004). The dashed blue, dot-dashed 
green, dotted red and solid pink curves correspond to the boundary condition 
m = —0.03, 0.03, —0.25, and 0.12 respectively. The dashed black line is the 
spinodal curve, on which the scalar curvature diverges. 



where a = a//i^ + m? /2T . The scalar curvature of eq.© for the metric of eq. ffTSj) . 
in the limit m — )■ (which is our region of interest), is given by i? = A/B, where 

A = h -2T {Ah^ + 4r + 1) Sinh (^^^ + 4 (/i^ + 2T) Tanh (^-^^ + 3hSech^ 

- 2T^ + 2h\4:T{T - 2) - 1) + 2T {4h\T + 1) + T) Cosh (^^^ 
B = 2^ (ranh (^^^ - 2h^ 

The scalar curvature diverges at Tanh^ = 2h, defining the phase boundary, a 
result that matches with that obtained in [13]. To understand the behavior of 
geodesies in this model, we set h = 0.2, which implies the critical temperature 
T = 0.236. Numerical solutions of this geodesic equations close to the critical 
point are plotted in fig.(|3]). Here, we have taken all the geodesies to start from 
(T, m) = (0.3,0.004). The dashed blue, dot-dashed green, dotted red and solid 
pink curves correspond to m = —0.03, 0.03, —0.25, and 0.12 respectively. Also 
shown in dashed black is the spinodal curve, i.e the locus of divergence of the 
scalar curvature arising out of the metric of eq. (fT8l) . We find that the geodesies 
show the same turning behavior as in the other mean field models discussed in 
the previous subsection. We also note that in the limit of T — )■ 0, gmm diverges 
and qtt 0. Numerical solutions seem to become somewhat unreliable in this 
limit, and we will not discuss them. 

Having elucidated the nature of geodesies in classical systems exhibiting phase 
transitions at non-zero temperatures, we finally move to quantum phase transi- 
tions at zero temperatures. 
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3 Geodesies in QPTs : The Transverse XY Spin 
Chain 



Information geometry of QPTs has been well studied of late, starting from the 
work of There are, however, very few systems to which this can be meaning- 
fully applied, since the definition of the geometry (from eq.(jl])) requires complete 
knowledge of the many body ground state, which may be difficult to obtain ex- 
cepting for a few exactly solvable system, like the transverse field XY spin chain. 
Even when such ground states are obtainable, as in the Dicke model of quantum 
optics, explicit calculations might be prohibitively difficult due to algebraic com- 
plications. We will base our calculations on the transverse XY model, for which 
the information metric was obtained in [4J. 

To recall, for the transverse XY spin chain, the Hamiltonian with {2N + 1) 
spins is 



H- 



XY 



j=-N 



(20) 



where the cr*, i = x,y,z are Pauli matrices, 7 is an anisotropy parameter, h is the 
magnetic field, and the Planck's constant has been set to unity. The information 
metric for this model has been calculated in [4] and in the thermodynamic limit, 
the line element, in the region < 1, 7 > (the ferromagnetic phase) is given 
by 

= I67 (1 - /^2) + I67 (1 + 7)' ^ ' 

QPTs occur on the lines 7 = 0, < 1 (the anisotropic transition line), and 
\h\ = 1, (the Ising transition lines), where the spectrum of the theory becomes 
gapless. Information geometry is however very different for these two transitions. 
Whereas the scalar curvature (calculated from eqs.([SD and fl?Il) ) diverges on the 
line 7 = 0, it is regular on the lines \h\ = ±1. For this model, the geodesic 
equations are 

^ + T^-T^ ' ^~ 27(1+7) +27(1-/.^)=° ^^^^ 

where, as before, the dot represents a derivative with respect to the affine param- 
eter r. Also, the normalization condition implies that 

167(1-/12)^167(1 + 7)2 ^ ' 

Before attempting to solve the coupled non-linear equations of eq. (l22|) . let us 
look at a special case. The first of eq. (!22|) is satisfied hj h = constant and hence 
constant h lines are geodesies. To find 7 as a function of the affine parameter 
in this case, we substitute h = in the second of eq. (|22l) and in eq.(l23l). Then 
it is seen that 7 = Tan^ (2 (r — tq)), where tq is a reference value for the affine 
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Figure 4: Numerical solution for a geodesic curve with (/i, 7) = (0.96,0.1) and 
(/i, 7) = (—0.0857, 1.35) on the /i — 7 plane. The geodesic is confined to a single 
phase region. 

parameter. Thus, for the constant h geodesies, 7 is always positive, i.e these 
geodesies do not cross the phase boundary at 7 = 0. Rather, they turn back on 
touching that line. This should be contrasted with the m = constant geodesies 
of the CW model, which, as we have said is not fully understood. 

To solve the equations in eq.( l22|) in general, we adopt a numerical procedure 
analogous to what we have done before. As an illustration, we solve for these 
equations with the initial conditions (/i, 7, h) = (0.96, 0.1, —0.0857). The solution, 
plotted on the h — j plane parametrically, with the affine parameter r, is shown 
in fig.dl]). Clearly, the geodesic is confined to a single phase, and does not cross 
the phase boundaries, as in CPTs. It is not difficult to check this analytically by 
expanding the metric near the lines 7 = and h = ±1. 

4 Conclusions and Discussions 

In this paper, we have studied four model systems that exhibit phase transitions, 
in the thermodynamic limit. The Van der Waals model, the Curie- Weiss mean 
field model of ferromagnetism and the infinite Ising ferromagnet exhibit CPTs at 
finite temperature. The transverse XY spin chain shows a QPT at zero tempera- 
ture. For all these models, we performed the most general analysis of geodesies in 
the parameter manifold. Such an analysis has not appeared in the literature be- 
fore. In the process, we have established the information geometry of the infinite 
Ising ferromagnet. We have solved the geodesic equations for all these models in 
full generality, by choosing a starting point (i.e coordinates) in the manifold, and 
imposing initial conditions on its derivatives with respect to the affine parameter. 
In this way, we are able to trace out the geodesies, and study their behavior near 
second order critical points. This complements and extends the results of [6], [7] 
in a non-trivial way. 

Our main conclusion here is that purely from a geometric perspective, geodesies 
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near critical points show universal behavior, although the physical nature of the 
phase transitions are widely different. We have also established that geodesies are 
confined to a single phase. We believe that these results are model independent, 
and should be true for any model of CPTs or QPTs. 

It might be interesting to study geodesies in the context of scaled equations 
of state for classical fluid systems, and also for some other models that exhibit 
QPTs. In particular, in the context of CPTs, it is an interesting question to 
ask if geodesies show any special behavior at or near the Widom line, which 
is a continuation of the co-existence curve, along which the correlation length 
maximizes. We leave such a study for the future. 
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